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abstract 

I review the recent progress in small x physics, concentrating on the topics 
relevant to the BFKL evolution. 



1 Introduction 

It is known that radiative corrections in perturbative QCD (PQCD) produce 
large logarithms, such as In Q in the kinematic region with large momen- 
tum transfer Q and a<jln(l/x) in the region with small Bjorken variable x. 
These large logarithms, spoling the expansion in the coupling constant as, 
must be organized in some way. To sum the various logarithms contained in 
a parton distribution function to all orders, four well-known evolution equa- 
tions have been proposed: The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equation sums InQ for intermediate x, the Balitskii-Fadin- 
Kuraev-Lipatov (BFKL) equation [0 sums ln(l/a;) for a small x, and the 
Ciafaloni-Catani-Fiorani-Marchesini (CCFM) equation 0, appropriate for 
both large and small x, unifies the above two equations. In the region with 
both large Q and small x many gluons are radiated in scattering processes 
with small spatial separation among them, and a new effect from the an- 

^talk presented at the International Symposium on QCD Corrections and New Physics, 
Hiroshima, Japan, 1997. 
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nihilation of two gluons into one gluon becomes important. The nonlinear 
Gribov-Levin-Ryskin (GLR) equation [Q takes into account this effect. All of 
the above equations can be applied to QCD processes at small x, though they 
give different evolutions in x of a parton distribution function. The DGLAP 
equation is applicable, since the relevant splitting functions show the desired 
X dependence. In this talk I will concentrate on the BFKL equation. 

2 The BFKL Evolution 

First, I list some important features of the BFKL evolution, which motivate 
the studies reported in the following sections. 

2.1 The BFKL Equation 

In the small x region gluon contributions dominate the structure functions 
involved in deep inelastic scattering (DIS) of a proton, such as F2(x, Q^), Q 
being the momentum transfer form the virtual photon. Since the longitudinal 
component ^p^ and the transverse component kx of the gluon momentum 
may be of the same order, the transverse degrees of freedom of the gluon 
must be taken into account, leading to the /cr-factorization theorem 0, 

Q') = f / ^H{x/t Q, kT)F{C, kr) , (1) 

Jx t, J TT 

where H is the hard photon-gluon scattering subamplitude, and F the unin- 
tegrated gluon distribution function. If x is large, the fcr dependence in H 
is negligible. The integration of F over kx then gives the gluon density G, 

^G{^,Q')= r^F{^,kT), (2) 

JO TT 

and Eq. (|T]) reduces to the conventional collinear factorization formula. 

It is known that leading logarithmic corrections asln(l/a;) are produced 
from a ladder diagram, where the rung gluons obey the strong rapidity or- 
dering, 

1/1 > Z/2 > ... > Vn , (3) 

with Hi the rapidity of the rung gluon closest to the proton, and i/n the rapid- 
ity of the rung gluon closest to the hard scattering. One further reggeizes the 
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ladder gluons by assuming that their invariant masses are given hj k"^ = —k^. 
Summing the Reggeon ladders to all orders, one obtains the BFKL equation 
that describes the evolution of F at small x [01, 



d 

x—f{x, kr) = ask^ 



f{x,k'j) - fix.kr) , fix.kr) 



I 



I 



+ 



{Ak'^-k^yi^ 



with /(x, kx) 
colors. 



(4) 

k'^F{x, kx) and = Ncas/n, Nc = 3 being the numebr of 



2.2 Solution to the BFKL Equation 

Using the Mellin transform 



fix,uj)= / ikY^-'f{x,k)de 
Jo 



1 rc+ioo _ 

f{x,k) = — {erf{x,u)du 

ZiTl Jc-ioo 



Eq. (^) in the uj space is written, for fixed a^, as 

d - 

where the trasformed BFKL kernel K is given by 

K{uj) = «,[2V^(1) - - - to)] , 

V'(w) = l-lnr(u;) , = -7^ , 

auj 

jE being the Euler constant. K{uj) has a maximum at u 
which it can be expanded as 

1 



1/2, 



with the constants A = a^41n2 and A" = as28C(3), C(3) = 1.202. 
tuting Eq. ( p!OD into (^, one solves for /, and the inverse transform 
gives 



f{x, k) 



xy^ /(xo,l/2)fc 
xo) [27rA"ln(xo/x)]V; 



exp 



(In A;^ — In A; 
2A"ln(xo/x) 



2\2 



(5) 
(6) 

(7) 

(8) 
(9) 
around 

(10) 

Substi- 
Eq. d) 

(11) 
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where 



In k"^ = i-^ 



InF {xo,^ + iu 



(12) 

u=0 



is the center of the Gaussian broadening in Infc^. xq is usually chosen as 0.1, 
below which the BFKL evolution begins. 

Equation (|TTp exhibits the following features: 

1. It possesses a power-law rise as with A ~ 0.5 for reasonable values 
of as, which is attributed to hard pomeron exchanges (a pomeron can be 
regarded as a color-singlet Reggeon ladder). Recent HERA data of the 
DIS structure function F2 confirm this feature. 

2. It does not depend on Q, implying that F2 is insensitive to the variation 
of the momentum trasfer. However, the rise of the data is characterized by 
A with a stronger Q dependence: A ~ 0.2 for ~ 4 GeV^ and A ~ 0.5 
for ~ 50 GeV^. At very low Q, A ~ 0.08 corresponds to soft pomeron 
exchanges. 

3. The transverse momentum diffuses into a nonperturbative region 



at small x. Since the Gaussian width ^2X" \n{xo/x) increases with 1/x, the 
distribution of the gluon in In fc^ certered at In k'^ becomes broader. The 
nonperturbative region with small k'^ then gives essentail contributions to 
F2, and PQCD is not reliable. 

4. The power-law rise of / renders F2 and the DIS cross section a violate 
the unitarity bound a < const, x ln^(l/x). 

The above features have stimulated intensive studies of the BFKL evolu- 
tion, whose progress will be reported below. 

3 How to explain data? 

To explain the data of the structure function F2, the DGLAP equation should 
be combined in some way, which introduces the Q dependence through the 
InQ summation. 

3.1 The DGLAP Equation 

One may employ the DGLAP equation for the gluon density G directly IQ, 
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The required rise of G at small x is the result of the splitting function 

1 



1 

7+1 



: 1 



which diverges at 2; ^ 0. It can be shown that G evolves according to 
xG{x, Q^) oc exp 









In — In — 




X to. 



(14) 



(15) 



with the variable t(o) = \n{Q^Q-j / Aqq^)) /3o the first coefficient of the QCD 
beta function. 

For Ql = 4: GeV2, Eq. (|T|) gives A ~ 0.11 for ^ 15 GeV^ and A ~ 0.15 
for ~ 30 GeV^. Obviously, these values of A are too small to explain the 
rise of F2. Hence, a steeper input G{x, Ql) must be adopted to compensate 
the slower DGLAP evolution. If choosing Ql = 0.3 GeV^, A ~ 0.4 for 
~ 20 GeV^ is large enough, and a flat input serves the purpose [§]. 
However, it has been criticized that the PQCD evolution around this low Qo 
is not reliable 0. 

The fcr factorization shown in Eq. (^ is in fact equivalent to the collinear 
(mass) factorization |^ , on which the DGLAP equation is based. Expressing 
the solution to Eq. (|^ in the moment space, one obtains 



/jv(u;) oc 



with the definition 



fN{uj) 



N -Kiuj) 



dxx^ ^f{x, iu) . 



(16) 



(17) 



Equation (|^) implies a pole of fN{^) at = 77V, which satisfies K{'j]y) = N, 
and the parametrization of -F/v(/c) from Eq. (0) p[. 



FNik) 



In 



7jv 



— G'iv(/i 



(18) 



H being an arbitrary factorization scale. 7iv is the BFKL anomalous dimen- 
sion, whose perturbative expansion is given by 



Inserting the above expression into Eq. (|1|) in the moment space, one arrives 
at 

- CM'/iJ.')G^if,') , (20) 



F2,n{Q 

with the coefficient function 



In 



(21) 



It can be shown that Gis[{ii^) is related to the iV-th moment of the gluon 
density G defined in the modified minimal subtraction scheme 0, and the 
BFKL anomalous dimension 77V is the same as the gluon anomalous dimen- 
sion 7c,c,,Ar, the A^-th moment of Pgg in Eq. (p!^, up to leading logarithms. 
Obviously, Eq. ([20|) is equivalent to the coUinear factorization formula for F2. 



3.2 Next-to-leading Logarithms 

As a more complete analysis, one includes the singlet quark contribution 
through the DGLAP equation 
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(22) 



where Q denotes the singlet quark distribution function. The anomalous di- 
mensions 7gg and 7,3, giving the next-to-leading-logarithm summation, leads 
to a steeper rise of G at small x compared to the DGLAP evolution in 



Eq. (|T5D, which takes into account only the leading gluon contribution. The 
analysis concludes that a flat input G{x, Ql) at = 4 GeV^ is preferred. 

The above conclusion has been justified in an alternative way. One ex- 
tracts the evolution equation for the gluon decaisy G from Eq. (|^) 



d ^ d 2 
+ 5— - 7 



(23) 



with 




(24) 
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where the variable t has appeared in Eq. (pISl). With the further variable 
changes 



a 



Eq. (|23D is easily solved to give 



(25) 



G{a, p) oc exp \ 2'^a — 5— 

Equation ( ]26| ) exhibits two asymptotic scaling laws: 



InG oc 0" 
InG oc const. 



for fixed p , 
for fixed a . 



(26) 



(27) 
(28) 



The plots of the data of lnF2 with respect to cr and p indeed confirm the 
double scaling behavior at large a and p. On the other hand, different in- 
puts of G{x, Qq) lead to different slopes of the scaling. It was found that 
the predictions from the flat input match the data better. Because of the 
above more complete analyses, one is tempted to conclude that soft pomeron 
exchanges dominate in the low Q region. It was also argued that the BFKL 
pomeron appears only for []TT] 



In ( - ) > 



asiQ) 



20 



(29) 



3.3 The CCFM Equation 

At last, one may resort to the CCFM equation, which embodies both the 
InQ and ln(l/x) summations. It is written as 

F{x,pT,Q)=F^'^\x,pT,Q)+ f dz [ -^e{Q-zq)AsiQ,zq) 



xP{z,q,pT)F{x/z,\^T+ {I- z)cilq) , (30) 



with the function 

P = ttsipr) 



h Ans{z, q,PT)- + z{l - z) 

1 — z z 



(31) 
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similar to the gluon splitting function Pgg in Eq. The so-called "Su- 

dakov" exponential and the "non-Sudakov" exponential A^s are given 
by 



^s{Q,zq) = exp 



—a 



(lp2 rl-pT/p dz' 



A 



NS[Z,q,PT) 



exp 



^0 dz' 



Pt 



dp^ 



Z' J{z'q)2 



(32) 



where the upper bound zq of the variable z' takes the values 0, |T^ 

1 if 1 < iPr/q) 

Zq = Pt/ q li z < {pr/ q) <l 
z if {pt/ q) <z . 



(33) 



A5 collects the contributions from the ladder diagrams with rung gluons 
obeying the strong angular ordering, and is the result of the In Q summation. 
Those gluons which do not obey the angular ordering are grouped into A^vs- 
In the small x region one adopts the approximate version of the CCFM 
equation |T^, 



F{x,pT,Q) = F^'\x,pT,Q) + as{pT) t- j ^e{Q - zq)e{q - f^) 

Jx Z J TXQ 



Txq^ 



xANs{z,q,PT)F{x/z, Ipt + (1 - z)q\,q) . (34) 



where the extra function 0{q — jj) introduces an infrared cutoff of the variable 
g, and F'-'^^ is a flat nonperturbative driving term. 

The CCFM evolution with x has been extracted from Eq. (0) given 

by 



F oc exp[2a,ln(g//i)ln(l/x)] . 



(35) 



It shows a steeper rise at small x than the DGLAP evolution in Eq. (|T^), 
though it also considers only the leading gluon contribution, and a more 
Q-dependent rise than the BFKL evolution in Eq. (ITTl). This is the reason 
the CCFM equation can explain the data well ||12[. The steeper CCFM rise 
compared to the DGLAP rise is attributed to the different choices of the 
argument of the running coupling constant in the splitting function, which is 
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Pt in the former (see Eq. (pT])) and Q in the latter. Therefore, asipr) does 
not run in fact, as the variable q is integrated over in the CCFM equation 
(|51[). While as{Q) runs, when solving the DGLAP equation. Because of 
<^s{Pt) > CisiQ), the CCFM evolution is stronger. 



4 Where is the hard pomeron? 

In the previous section it has been stated that soft pomeron contributions 
dominate in the low-Q region, and the high-Q behavior of the DIS structure 
function is the consequence of the DGLAP evolution plus a nonperturbative 
input. Though the BFKL equation gives an equally good explanation of the 
data at large Q, it is still worthwhile to ask where the clear evidence for the 
BFKL (hard) pomeron is. To answer this question, many proceses have been 
proposed. 

4.1 Two-jet Processes 



Consider the two-jet production from hadron-hadron collisions [14 



p + p^2 jets + X . (36) 

If tagging the jet with longitudinal momentum fraction Xi, transverse momen- 
tum piT and rapidity t/i, i = 1, 2, the momentum fractions of the initial-state 
partons will be fixed at Xi due to the kinematic relation 

= (Pit/Vs) exp(?/j) , (37) 

s being the center-of-mass energy, when the rapidity gap Ay = yi—y2 is large. 
It was argued that many minijets are radiated from the t-channel exchanged 
gluon between the initial-state partons at large Ay, such that the cross section 
(P a / {dyidy2) increases with Ay. This rise is definitely attributed to the 
BFKL pomeron contribution, since the nonperturbative parton distribution 
functions take the values at x = Xi (no evolution in x) after the jet tagging. 
However, an explicit evaluation of (Pa / {dyidy2) showed that it, contrary to 
the expectation, decreases with Ay in the available range of s [|15||. The 



reason is that Xi approaches unity for finite piT and s as Ay increases, and 
thus the parton distribution functions vanish. 
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An alternative was then suggested, where the angle decorrelation between 



the two tagged jets is measured |jT5|. The radiation of the minijets washes 
out the correlation of the two jets with angle = vr at a rate that increases 
with Ay. The explicit analyses of the decorrelation have confirmed this 
expectation . 



4.2 One Jet in DIS 



A similar process with one forward jet in DIS was proposed |T^. By tag- 
ging the jet with momentum fraction Xj and transverse momentum kjT-, the 
momentum fraction of the initial-state parton is fixed at Xj. kjT is chosen 
as being large enough to avoid the kr diffusion into the nonperturbative re- 
gion. The t-channel exchanged gluon radiated from the parton emitts many 
minijets, such that the structure function F2 behaves as 



dFoix) 



X 



d\n{l/xj)dkjT 



oc as{kjT)xjf{xj) — 



•^3 , 



(38) 



Note that the parton distribution function 



/ 



Q) 



(39) 



is evaluated at a; = Xj, and the nonperturbative ambiguity is removed. The 
rise x~^ is then clearly identified as the consequence of the BFKL evolution. 
Though it is challenging to tag a forward jet, experimental studies have been 



attempted with encouraging results |T7 



Instead of tagging a forward jet, it has been suggested to tag a forward 
photon [0. The advantage of this process is that it is cleaner. However, it is 
suppressed by the coupling constant a em compared to the previous process, 
and the photon needs to be isolated from vr''. 



4.3 Levin's Viewpoint 

Though all the above analyses and experiments have not yet fully confirmed 
the existence of the BFKL pomeron. Levin argued that its evidence had been 



contained in the current data 19|. To demonstrate his viewpoint, an average 
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of the BFKL anomalous dimension (7) and the strength of the shadowing 
corrections k are defined by 

The BFKL anomalous dimension and the moment of the gluon density 
Gn that obeys the DGLAP evolution have been introduced in Sect. 3.1. 
(7) can be regarded as measuring the BFKL content in the DGLAP gluon 
density employed in the explanation of the DIS data |^. If (7) is located 
between 1/2 and 1, it is claimed that the gluon density contains significant 
BFKL pomeron contribution. From the GLR equation (see Eq. (^Tj) below), 
it is easy to observe that k corresponds to the nonlinear term, describing the 
effect of two gluon annihilation into one gluon. R is the correlation length 
of the radiated gluons. For k > 1, the shadowing correction is considered to 
be large. 

A simple analysis showed that the BFKL region between (7) = 1/2 and 
(7) = 1 in the ln(l/x)-ln plane indeed penetrates the region which has 
been explored by the HERA experiments. However, the BFKL region is 
completely located in the region with k > 1. Hence, Levin concluded that 
the BFKL pomeron contribution has been contained in the current HERA 
data, but is suppressed by the shadowing correction, and thus unobservable. 

5 How to Recover Unitarity? 

As stated in Sect. 2, the power-law rise violates the unitarity bound, 
implying that the BFKL equation derived from the leading-twist and leading- 
logarithm approximation needs to be corrected. The attempts to include 
higher-twist and next-to-leading logarithmic contributions have been made. 

5.1 The GLR Equation 

In the region with both large Q and small x, many gluons are radiated 
by partons with small spatial separation among them. A new effect from 
the annihilation of two gluons into one gluon is then essential. Taking into 
account this effect, the BFKL equation is modified by a nonlinear term. 
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leading to the GLR equation 



d'^xG{x,Q'^) _ , 2\ r ryf /n2M2 



The constant 7 = 81/16 is regarded as the effective couphng of the annihi- 
lation process, and the radius R characterizes the correlation length of the 
radiative gluons as mentioned before. It is obvious that the second term is of 
higher twist, but becomes important as R is small. The minus sign in front 
of it implies that the annihilation decreases the number of gluons, and that 
the rise of the gluon density at small x due to the first term might saturate. 

This equation is, however, nonlinear, and the values of R depend on 
models. It is set to 2 GeV~^ for the model with the radiative gluons concen- 
trated at hot spots, and 5 GeV~^ for the model with the gluons uniformly 
distributed in a hadron. 

5.2 Next-to-leading Logarithms 

To recover the unitarity in the framework of the BFKL equation, it was 
proposed to compute the kernel that includes next-to-leading ln(l/a;), which 
arise from relaxing the strong rapidity ordering of rung gluons. That is, the 
contribution from the region with, for example, 

1/1 -1/2 >Z/3 > - >l/n , (42) 

should be computed. The virtual and real corrections have been obtained in 
21] and in [^, respectively. However, an explicit expression of the BFKL 



kernel including next-to-leading logarithms has not been extracted from the 
above analyses (the phase-space integrals for the amplitudes of real gluon 
emissions given in were not performed). Hence, no concrete conclusion 
on the unitarity has been drawn. Since the involved calculations are very 



tedious, I will not go into the details. Readers are referred to Ref. ||2^ for 
summary of the progress on this subject. 

5.3 Multiple Pomeron Exchanges 

Multiple pomeron exchanges, as higher-twist contributions, have been stud- 



ied by means of colored-dipole scatterings [^]. Consider an onium, which is 
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the QQ bound state with Q a heavy quark. The gluon emission inside the 
onium, from the viewpoint of color flow, can be regarded as being composed 
of the Q-Q pair in the color octet state. Many gluons then imply many quark 
pairs. The emitted quarks and valence quarks, combined among each other, 
form many colored-dipoles in the oniun. In this picture, the onium-onium 
forward scattering is then formulated as multiple dipole scatterings. One 
of the advantages of this process is that the onium size provides a natural 
infrared cutoff, and thus resolves the problem of the kx diffusion. 

To compute the cross section, one defines a dipole density n{x, Y) for 
an onium, with x the transverse size of a dipole and Y the rapidity. At 
high energy, it has been shown that n increases with Y due to the same 
mechanism as of the rise of the gluon distribution function with 1/x. That 
is, the equation that governs the evolution of n with Y is equivalent to the 
equation that governs the evolution of F with ln(l/x) However, the 

derivation of the kernel for the former is much simpler than that for the 
latter. Though the whole evolution kernels for n and for F are equal, the 
parts from real gluon emissions and from virtual gluon emissions differ. To 
obtain the kernel for F, all the higher-order diagrams are considered. For the 
kernel of n, only the diagrams with radiative gluons emitted before the hard 
onium-onium scattering are considered, since they are the diagrams that 
modify the dipole density for an initial-state onium. With this criterion, 
the number of diagrams needed to be computed is greatly reduced. This is 
another advantage of the above process. 

Since the formalism for dipole scatterings is simpler, it can be extended 
to include multiple pomeron exchanges. Recall that the t-channel gluon 
responsible for the hard scattering is reggeized at large rapidity gap, and two 
reggeons form a pomeron. The cross sections including single and double 



pomeron exchanges have been obtained, which exhibit a rise with Y ||25|. As 
to the unitarity, for which infinite many pomeron exchanges may be essential, 
no conclusive progress has been made following this vein. 

5.4 Reggeon Compound States 

The higher-twist contributions from multiple pomeron (reggeon) exchanges 
can be taken into account in another interesting approach . As explained 



in Sect. 2, a reggeized gluon is a two-dimensional object, since its prop- 
agator depends only on the transverse momentum. Hence, in high-energy 



13 



scatterings (3+1) QCD is reduced to (2+1) Regge theory |2^. Consider the 



compound state of reggeons distributed randomly in a x-y plane. It has 
been shown that only the nearest-neighbor interaction survives in the large 
Nc limit. Under this approximation, these reggeons can be arrayed into a one- 
dimensional space in the sequence where the interaction occurs between any 
two neighboring reggeons. Certainly, there exist two such one- dimensional 
spaces, characterized hy z = x + iy and z = x — iy. Then the (2+1) Regge 
theory is further reduced to a (1 + 1) ® (1 + 1) Schrodinger equation. It 
is trivial to observe that each of the one- dimensional systems is an exactly 
solvable Heisenberg model with spin s = 0. 

The Hamiltonian and the corresponding Schrodinger equation are written 

as 



Hn{z, z) = as [Hn{z) + Hn{z)] , (43) 
HnXn = EnXn , (44) 

where xn is the wave function of the A^-reggeon compound state with the 
eigenenergy E^. The eigenenergy can be identified as the intercept of the 
Regge trajectory. Therefore, by computing the eigenenergy of the compound 
state, one understands how the power-law rise of the structure function F2 
varies with the reggeon number A^: 

Hx,Q2) ^a;-^^ . (45) 

E2 = Q;s41n2 for the N = 2 case, which corrsponds to the BFKL pomeron 
exchange, was found. The tendency E2 > E^ > ... has been observed, in- 
dicating that the rise is indeed softened gradually as increases, though 
it is still power-like. Since the N ^ 00 case remains unsolved, a concrete 
conclusion on the unitarity is still not available. 



6 Summary 

In this talk I reviewed the recent progress in small x physics. Obviously, 
many issues need to be studied theoretically and experimentally with more 
effort. 
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